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Abstract 

Background: Tlie dynamic re-emergence of visceral leishmaniasis (VL) in south Europe and the northward shift to 
Leishmanio-free European countries are well-documented. However, the epidemiology of VL due to Leishmonio 
infontum in southeastern (SE) Europe and the Balkans is inadequately examined. Herein, we aim to re-evaluate and 
compare the population structure of L infantum in SE and southwestern (SW) Europe. 

Methods: Leishmania strains collected from humans and canines in Turkey, Cyprus, Bulgaria, Greece, Albania and 
Croatia, were characterized by the K26-PCR assay and multilocus enzyme electrophoresis (MLEE). Genetic diversity 
was assessed by multilocus microsatellite typing (MLMT) and MLM Types were analyzed by model- and 
distance- based algorithms to infer the population structure of 128L infantum strains. 

Results: L infantum MON-l was found predominant in SE Europe, whilst 16.8% of strains were MON-QS. Distinct 
genetic populations revealed clear differentiation between SE and SW European strains. Interestingly, Cypriot canine 
isolates were genetically isolated and formed a monophyletic group, suggesting the constitution of a clonal MON-l 
population circulating among dogs. In contrast, two highly heterogeneous populations enclosed all MON-l and 
MON-98 strains from the other SE European countries. Structure sub-clustering, phylogenetic and Splitstree analysis 
also revealed two distinct Croatian subpopulations. A mosaic of evolutionary effects resulted in consecutive 
sub-structuring, which indicated substantial differentiation and gene flow among strains of both zymodemes. 

Conclusions: This is the first population genetic study of L infantum in SE Europe and the Balkans. Our findings 
demonstrate the differentiation between SE and SW European strains; revealing the partition of Croatian strains 
between these populations and the genetic isolation of Cypriot strains. This mirrors the geographic position of 
Croatia located in central Europe and the natural isolation of the island of Cyprus. We have analysed the largest 
number of MON-98 strains so far. Our results indicate extensive gene flow, recombination and no differentiation 
between MON-l and MON-98 zymodemes. No correlation either to host specificity or place and year of strain 
isolation was identified. Our findings may be associated with intensive host migration and common 
eco-epidemiological characteristics in these countries and give valuable insight into the dynamics of VL. 




Parasites 
SiVectors 



* Correspondence: ksoteriadou@pasteur.gr 

laboratory of Molecular Parasitology, Hellenic Pasteur Institute, 127 Vasilissis 
Sofias Avenus, 1 1521, Athens, Greece 

Full list of author information is available at the end of the article 

O© 2013 Gouzelou et a!.; licensee BioMed Central Ltd. This is an open access article distributed under tine terms of the Creative 
BIoIVIGCI CGntrsI commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Gouzelou et al. Parasites & Vectors 2013, 6:342 
http://www.parasitesandvectors.conn/content/6/1/342 



Page 2 of 1 9 



Background 

Leishmania infantum is the main causative agent of hu- 
man visceral leishmaniasis (HVL) and canine visceral 
leishmaniasis (CanL) and also responsible for sporadic 
cases of cutaneous leishmaniasis (CL) across the Mediter- 
ranean basin [1,2]. The concurrent species Leishmania 
tropica extending rarely to Europe via Greece [1,3-5] and 
Leishmania major present in North Africa and Middle 
East, are also accountable for anthroponotic and zoonotic 
CL, correspondingly [1]. 

Based on multilocus enzyme electrophoresis (MLEE), 
the current reference method for characterizing and 
classif)^ing Leishmania strains [6], L, infantum MON-1 is 
the predominant zymodeme in all Mediterranean coun- 
tries [7-9], while the genetically close L, infantum zymo- 
deme MON-98 has occasionally been reported in Turkey 
[10], Cyprus [11], Egypt [12,13], Greece [5,14,15] and 
Portugal [16]. 

Dogs constitute the main reservoir host of L, infantum 
MON-1 [9] and probably MON-98 [5], while cats could 
be a secondary domestic L, infantum MON-1 reservoir 
[17,18]. Feline leishmaniasis is reported in countries highly 
endemic for CanL [19] and cats were recently found 
capable of transmitting MON-1 parasites to a proven 
L, infantum vector, Phlebotomus perniciosus [17]. The 
large number of sub-clinically infected dogs is a major 
veterinary and public health problem in south Europe 
[1]. Although clinically affected dogs are more prone to 
infect sand fly vectors, those being sub-clinically infected 
and seronegative can also transmit Leishmania parasites 
to sand flies and, therefore, contribute to the parasites' 
maintenance [20-24]. 

Further to this, the distribution range of CanL has 
surpassed the limits of south Europe and is spreading to- 
ward northern European countries, presently considered 
non-endemic. In Germany, there is an estimate of 20,000 
infected dogs that were either imported from the Mediter- 
ranean basin or infected after travelling to this region 
[25,26]. Furthermore, the presence of a L, infantum vector 
(P. perniciosus) in Germany and Switzerland has been re- 
cently reported [27]. The northward spread of CanL is 
also apparent in northern Spain, where new foci of CanL 
and the presence of its vectors {Phlebotomus ariasi and 
P, perniciosus) have been reported. In Italy, canine inci- 
dence has reached 30.3% in West Liguria [28], and in 
France seroprevalence fluctuates between 3.2% in Alpes- 
Maritimes and 26.5% in Corsica [29]. 

Concurrently, a considerable increase in VL cases due 
to L, infantum MON-1 is documented in traditionally 
endemic southeastern (SE) European countries. Since 
the first reported HVL cases in the islands of Spetses 
(1879) and Hydra (1881) [30], leishmaniasis is endemic 
in most islands and coastal regions of Greece and consti- 
tutes a significant veterinary and public health concern. 



During 1994-2001, HVL has re-emerged in inland re- 
gions of Greece, where seroprevalence reached 33.4% in 
Attica and 12.6% in Epirus [31]. During the same 8-year 
survey, the highest CanL seroprevalence in clinically af- 
fected dogs was found in Attica (48.4%) followed by Epirus 
(45.4%) and central Macedonia (40%) [31]. Of additional 
concern are the high seropositivity rates found in clinically 
healthy dogs across the mainland of Greece. According to 
recent studies, seroprevalence in healthy dog populations 
ranged between 22.4-30.1% in Attica, 21-24.4% in Epirus, 
12.3%-17% in central Greece and 10.8%-20.5% in central 
Macedonia [31,32]. Among all Greek islands most dra- 
matic changes have been witnessed in Crete, where CanL 
increased from 0.3% to 19.8% during 1990-2009 and HVL 
prevalence from 0 to 1.33 per 100,000 inhabitants [4,5]. 

In south Cyprus, HVL and CanL do co-exist, but 
show different epidemiological patterns. In contrast to 
other endemic Mediterranean countries, in Cyprus two 
separate concurrent disease patterns involve CanL due to 
L, infantum MON-1 [11,33] and HVL due to Leishmania 
donovani MON-37 [34-36]. CanL was widespread by 
1935 and constituted a major veterinary health problem 
until 1945 [33]. Since the completion of subsequent eradi- 
cation campaigns against malaria and echinococcosis, 
CanL gradually re-emerged in coastal areas and is now 
widespread and highly endemic. In particular, the overall 
canine incidence has risen from 1.7% to 14.9%- 19.6% 
[11,33]. Currently, the highest prevalence rates found in 
all prefectures of south Cyprus are 33.3% (Larnaca), 32.4% 
(Limassol), 30.8% (Paphos) and 20% (Nicosia), while 24.3% 
of the PCR positive and 9.2% of the seropositive dogs are 
subclinically infected [11]. HVL due to L, infantum MON- 
1 has been absent since 1987 [33] and there are no official 
reports of infected natives even in areas highly endemic 
for CanL in south Cyprus. Yet, an increased incidence of 
human leishmaniasis possibly due to L, infantum has been 
reported in the northern part of the island [11,37]. 

In Turkey, the presence of leishmaniasis was first doc- 
umented in the early 19th century in Trabzon (eastern 
part of Black Sea region) [38]. Nowadays, CanL is wide- 
spread across the country (0.7%-33.3%) [10,39] and highly 
endemic along the Aegean and Mediterranean coasts 
(11-30.4%) [10]. Childhood VL is also well-established, 
rising in recent years [10] and currently highly endemic in 
Central Anatolia, the Aegean and Mediterranean Regions 
[38]. However, the most recent data from Denizli province 
(Aegean Region) showed that all children were seronega- 
tive whereas adult HVL had increased significantly [10]. 
As in neighboring Cyprus, parasites belonging to a genet- 
ically distinct L, infantum non MON-1 group have been 
reported to cause VL or CL in Turkey [34,36]. 

In Croatia, HVL and CanL existed across the Dalmatian 
coasts (Split, Makarska, and Dubrovnik areas) until 1957. 
The majority of reported HVL cases (90.4%) involved 
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children under 10 years old and 4% of these cases were 
fatal [40]. The disease occurred only sporadically for 
30 years due to mass spraying with insecticides for malaria 
control, but dramatically re-emerged during the war and 
postwar periods. During the past two decades, human in- 
cidence has reached 0.4/10^ of the population in southern 
Dalmatia [40] and is still considered a paediatric disease as 
64% of cases involve children [41]. Reports on HVL cases 
suggest that existing differences regarding the clinical 
course of infection and/or disease outcome may be associ- 
ated, amongst other factors, with variations of parasite 
species [42]. CanL is endemic in the city of Split (14.7%) 
and across the enzootic Split-Dalmatia County (7.1%- 
42.8%), while 45.7% of seropositive dogs were found to be 
sub-clinically infected [43]. 

In Albania, the incidence and geographical spread of 
HVL and CanL has gradually increased since World 
War II. However, the disease sharply re-emerged during 
1997-2001 with an average 173 HVL cases/year and 2.8/ 
10^ population morbidity, which exceeded 1/10^ lethality 
in active foci across north and central Albania [44]. En- 
demic zones are no longer restricted to northern regions 
or coastal and hilly areas. Children comprise 70-87.7% 
of HVL cases with an annual incidence of 25/10^ (0-6 year 
age group) [44-46]. Adult HVL cases and WW I Leishmania 
co-infections are also recorded [45]. Since 1977, CanL is 
widespread across the country with 16-17% seroprevalence 
recorded in Tirana and several other districts [44,46]. 

In Bulgaria, canine seroprevalence has gradually de- 
creased from 81% in 1941 to 13.8% in 1993 [47,48]. How- 
ever, during 1988-2004 an HVL outbreak occurred in 
south Bulgaria and an unexpectedly high morbidity rate 
(50.7%) was noticed in Stara Zagora. In 2005, two autoch- 
thonous HVL cases were reported in the same region, but 
all stray dogs tested were found seronegative [49]. Con- 
versely, in 2006 two breed dogs, which were not imported 
and had never travelled outside Bulgaria, were diagnosed 
with CanL [47]. L infantum is incriminated for human 
and canine disease [49], underreporting is considered mild 
and the estimated annual VL incidence ranged between 8 
and 12 cases/year during 2004-2008 [50]. 

The epidemiology of VL in these SE European and 
Balkan countries has been either inadequately or never 
examined so far. Multilocus microsatellite typing (MLMT), 
a powerful tool for phylogenetic and population genetic 
studies in Leishmania [51-55], previously showed that 
L infantum MON-1 strains from Greece differentiate from 
those of Western Mediterranean [15] and North Africa 
[56,57]. In addition, there is no indication of subdivision 
between Turkish and Greek strains or within Greek strains 
as well as no differentiation between MON-1 and MON- 
98 strains from Greece [15]. These findings prompted us 
to initiate a comprehensive population genetic study in SE 
Europe. For this, 109 Leishmania strains (25 human and 



84 canine isolates) were collected from Turkey, Cyprus, 
Bulgaria, Greece, Albania and Croatia, and characterized 
by MLEE and K26 typing. The analysis of L infantum 
strains by MLEE has demonstrated a broad enzymatic 
polymorphism within this species but has limited epi- 
demiological applicability in Europe due to the predomin- 
ance of MON-1 [6]. The /<26-PCR assay, which is specific 
for the L donovani complex, allows simultaneous discrim- 
ination between L infantum MON-1 and more than nine 
L infantum non MON-1 zymodemes that are known to 
exist in Europe [14,58,59]. Notably, the 626 bp/<26 ampli- 
con corresponds only to zymodeme MON-1, while the 
940 bp K26 amplicon may correspond either to zymodeme 
MON-1 or MON-98. For microsatellite typing, a standard 
set of highly discriminative markers [52] was employed for 
all strains from SE and their MLMT profiles were com- 
pared to 19 previously analyzed strains of L, infantum 
MON-1 from southwestern (SW) European countries 
(France, Spain and Portugal) aiming to re-assess the 
genetic diversity and population structure of L, infan- 
tum in SE Europe. 

Methods 

Leishmania strains and study area 

Table 1 summarizes the characteristics of the 128 L. infan- 
tum strains included in the study. Fifty-four strains were 
isolated from 20 human and 34 canine VL cases in diverse 
geographic regions of the mainland and islands of Greece. 
Sixteen human and canine isolates were obtained from 
Turkey, mainly from the Aegean region. Of these, nine 
MON-1 or MON-98 strains from Greece and two Turkish 
strains of unknown zymodeme type (strains EPS and EP16) 
were previously analyzed in different studies [15,52]. 
Twenty-three canine isolates were collected from five pre- 
fectures of southern Cyprus and sixteen strains were iso- 
lated in Bulgaria, Albania and Croatia, mostly from canines. 
Nineteen previously analyzed L infantum MON-1 strains 
from France, Spain and Portugal [15,52] were incorporated 
for comparison. 

Parasite cultures, cloning and DNA extraction 

Leishmania promastigotes were cultured at 25°C in 
RPMI 1640 or Ml 99 cell growth medium supplemented 
with 10-20% heat-inactivated FBS, 5U/ml penicillin, 5 (ig/ 
ml streptomycin and 10 mM HEPES buffer. The number 
of in vitro passages was limited to a minimum to avoid 
contamination and genetic drift due to serial sub-culturing. 
Clonal populations from the mixed L. infantum MON-l/L 
donovani MON-37 culture of the CD44 strain were ob- 
tained as described elsewhere [36]. Leishmania DNA was 
extracted from parasite cultures using a classical phenol/ 
chloroform protocol [60] or the standard protocol for cul- 
tured cells using the NucleoSpin® Tissue kit (Macherey- 
Nagel). Extracted DNA was respectively dissolved in dH20 
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Table 1 Characteristics of the MS Leishmania infantum strains studied 



Origin 


Strains 


WHO code 


Periphery/Region- 
Prefecture/Province (District) 


Pathology^ 


mon'' 


K26 

ampiicon'^ 


MLMT 
profile*^ 


Structure 

K = 4'^ 


TR(16) 


EP3* 


MHOMm/1994/EP3 


Aegean-Izmir 


VL 


n.d. 


n.d. 


35 


pop2 




EP16* 


MCANm/1996/EP16 


Aegean-Manisa 


CanL 


n.d. 


n.d. 


36 


pop2 




EP43 


MCAN^R/2000/EP43 


Aegean-Aydin(Ku§adasi) 


CanL 


98 


870/940 


37 


pop4 (sub4A) 




EP44 


MCAN^R/2000/EP44 


Aegean-Aydin(Ku§adasi) 


CanL 


98 


870/940 


37 


pop4 (sub4A) 




EP49 


MCAN^R/2000/EP49 


Aegean-Aydin(Ku§adasi) 


CanL 


98 


870/940 


37 


pop4 (sub4A) 




EP50 


MCANm/2000/EP50 


Aegean-Aydin(Ku§adasi) 


CanL 


98 


870/940 


38 


pop4 (sub4A) 




EP53 


MCANm/2000/EP53 


Aegean-Aydin(Ku§adasi) 


CanL 


98 


870/940 


39 


pop4 (sub4A) 




EP55 


MCANm/2000/EP55 


Aegean-Aydin(Ku§adasi) 


CanL 


98 


870/940 


40 


pop4 (sub4A) 




EP56 


MCANm/2000/EP56 


Marmara-Bilecik 


CanL 


1 


870/940 


25 


pop2 




EP57 


MCAN^R/2001/EP57 


Aegean-Mugia 


CanL 


n.d. 


870/940 


25 


pop2 




EP58 


MCAN/TR/2001/EP58 


Aegean-Mugia 


CanL 


n.d. 


940 


41 


pop2 




EP68 


MCANm/2001/EP68 


Aegean-Aydin(Ku§adasi) 


CanL 


n.d. 


940 


41 


pop2 




EP88 


MCANm/2003/EP88 


Black Sea-Corum 


CanL 


1 


626 


42 


pop3 (sub3B^) 




EP89 


MHOMm/2003/EP89 


Aegean-Manisa(Ala§ehir) 


VL 


n.d. 


870/940 


43 


pop2 




EP126 


MCANm/2005/EP126 


Aegean -Mugla(Bodrum) 


CanL 


1 


870/940 


13 


pop2 




EP134 


MCANm/2007/EP134 


Mediterranean-Hatay(Antakya) 


CanL 


1 


626 


44 


pop4 (sub4C^) 


CY(23) 


CDl 


MCAN/CY/2001/CDl 


Limassol 


CanL 


1 


626 


29 


popl (sublA) 




CD13 


MCAN/CY/2003/CD13 


Limassol 


CanL 


1 


626 


29 


popl (sublA) 




CD21 


MCAN/CY/2004/CD21 


Paphos 


CanL 


1 


626 


29 


popl (sublA) 




CD22 


MCAN/CY/2004/CD22 


Limassol 


CanL 


1 


626 


29 


popl (sublA) 




CD34 


MCAN/CY/2005/CD34 


Limassol 


CanL 


1 


626 


30 


popl (sublB) 




CD37 


MCAN/CY/2005/CD37 


Paphos 


CanL 


1 


626 


29 


popl (sublA) 




CD40 


MCAN/CY/2005/CD40 


Paphos 


CanL 


1 


626 


29 


popl (sublA) 




CD43 


MCAN/CY/2005/CD43 


Paphos 


CanL 


1 


626 


29 


popl (sublA) 




CD44 cl.5^ 


MCAN/CY/2005/CD44CI.5 


Paphos 


CanL 


1 


626 


31 


popl (sublB) 




CD48 


MCAN/CY/2005/CD48 


Limassol 


CanL 


1 


626 


29 


popl (sublA) 




CD51 


MCAN/CY/2005/CD51 


Paphos 


CanL 


1 


626 


29 


popl (sublA) 




CD57 


MCAN/CY/2005/CD57 


Nicosia 


CanL 


1 


626 


29 


popl (sublA) 




CD62 


MCAN/CY/2005/CD62 


Paphos 


CanL 


1 


626 


29 


popl (sublA) 




CD68 


MCAN/CY/2005/CD68 


Paphos 


CanL 


1 


626 


29 


popl (sublA) 




CD74 


MCAN/CY/2006/CD74 


Famagusta 


CanL 


1 


626 


29 


popl (sublA) 




CD63 


MCAN/CY/2005/CD63 


Limassol 


CanL 


1 


626 


29 


popl (sublA) 




CD66 


MCAN/CY/2005/CD66 


Nicosia 


CanL 


1 


626 


29 


popl (sublA) 




CD85 


MCAN/CY/2006/CD85 


Nicosia 


CanL 


1 


626 


32 


popl (sublB) 




CD87 


MCAN/CY/2006/CD87 


Paphos 


CanL 


1 


626 


29 


popl (sublA) 




CD88 


MCAN/CY/2006/CD88 


Paphos 


CanL 




626 


33 


popl (sublB) 




CD90 


MCAN/CY/2006/CD90 


Paphos 


CanL 




626 


29 


popl (sublA) 




CD91 


MCAN/CY/2006/CD91 


Paphos 


CanL 




626 


33 


popl (sublB) 




CDl 08 


MCAN/CY/2006/CD108 


Larnaca 


CanL 




626 


34 


popl (sublB) 


BG(1) 


BGl 


MCAN/BG/2008/BGl 


Blagoevgrad 


CanL 




940 


// 


pop2 


GR(54) 


GH2* 


MHOM/GR/2001/GH2 


Attica-Athens 


VL 




870/940 


1 


pop2 




GH5* 


MHOM/GR/2001/GH5 


Crete-Lasithi/Ag.Nikolaos 


VL 




870/940 


2 


pop2 




GH8* 


MHOM/GR/2001/GH8 


Attica-Athens 


VL 




870/940 


3 


pop2 



Gouzelou et al. Parasites & Vectors 2013, 6:342 
http://www.parasitesandvectors.conn/content/6/1/342 



Page 5 of 1 9 



Table 1 Characteristics of the MS Leishmania infantum strains studied (Continued) 



uHy 


)VlnU)VI/UK/zUU l/uHy 


Attica-Athens 


VL 


1 Q"7n/n/in 
1 o/U/y4U 


4 


pop4 (sub4C^) 


UH 1 U 


IVIHUIVI/LiK/zUU 1 /UH 1 U 


Attica-Athens 


\/i 
VL 


\ o/U/y4(J 


2 


pop2 


Un 1 1 


IVInUIVI/LDK/zUU l/Un 1 1 


Attica-Atnens 


\/i 
VL 


1 Q7n/Q/in 
1 o/U/y4U 


5 
J 


pop2 


UH 1 o 


IVInUIVI/UK/zUUo/uH 1 d 


Crete-Heral<lion 


VL 


yo o/U/y4U 


5 


pop4 (sub4A) 


run 


IVInUIVI/UK/zUU4/Unz 1 


Crete-Hera !<! ion 


\/i 
VL 


1 Q7n/Q/in 
1 o/U/y4U 


0 


pop2 


r\-\ll 


IVInwlVI/LDr\/ZUU4/Ur1ZZ 


Crete-Rethyinno 


\/l 
VL 


QQ Q7n/Q/in 
yo o/u/y4U 


£7 
J 


pop4 (^SUD4Aj 




IVInUIVI/Un/zUU4/Linzo 


Crete-Rethymno 


VL 


QQ Q7n/Q/in 

yo o/u/y4u 


7 
/ 


pop4 (sub4C^) 


Unz4 


KAur^KA/rD /onn/i /r^un/i 
)VlnU)VI/UK/zUU4/Unz4 


Attica-Atnens 


\/l 
VL 


1 DZD 


Q 
O 


pop4 (sub4C^) 


Unzo 


)VlnU)VI/UK/zUU4/UnzD 


Crete-Lasithi/Sitia 


\/l 
VL 


1 o/u/y4u 


Q 

y 


pop2 


Unz/ 


K /i 1-1 Pi N /I /r~D /inn /r"i-n7 
IVInUIVI/Un/zUUj/Linz/ 


Crete-Hera l<l ion 


\/l 
VL 


1 y4u 


1 n 
1 U 


pop2 


LLIVIzzoy 


IVIHUIVI/LiK/ 1 yyU/LA 1 Do/ 


Attica-Athens 


\/i 
VL 


1 o/(J/y4(J 


/ / 


pop2 


D A D A nr^ 
r Ar AUU 


)VlnU)VI/UK/ 1 yoU/r ArAUU 


Attica-Atnens 


\/i 
VL 


1 o/U/y4U 


5 
J 


pop2 




IVInUIVI/UK/zUUo/un4y 


North Aegean Sea-Lesbos/Mytilene 


\/i 
VL 


1 o/U/y4U 


Z 


pop2 


uHj I 


IVlnUIVI/bK/zUUo/bnj 1 


Central Greece-Euboea/Hall<ida 


\ /I 
VL 


1 940 


1 2 


pop2 


r"Uc:/l 


IVInUIVI/uK/zUUy/unj4 


Central Greece-Phthiotis(Lamia) 


\/i 
VL 


1 o/U/y4U 


7 3 
ID 


pop2 




IVInUIVI/UK/zUUy/uHjj 


Crete-Heraklion 


\/i 
VL 


1 Q7n/n/in 
1 o/U/y4U 


1 4 


pop2 




IVInUIVI/bK/zUUy/uHjo 


Peloponnese-Arcadia (Tripoli) 


\ /I 
VL 


1 o/(J/y4(J 


1 5 


pop2 




)VILAI\l/Un/zUU l/oUo 


Crete-Heraklion 


LanL 


QQ Q7n/Q/in 

yo o/U/y4U 


\ 0 


pop4 (sub4C^) 


uU4 


IVILAI\l/bn/zUU I/UU4 


Crete-Heraklion 


LanL 


y© o/U/y4U 


1 7 


pop4 (sub4C^) 




IVlLAN/uK/zUUz/LiUy 


Crete-Heraklion 


LanL 


\ o/U/y4(J 


1 8 


pop4 (sub4C^) 


UU 1 z 


IVILAIn/UK/zUUj/UU I z 


Crete-Lasithi/Sitia 


LanL 


1 Q7n/Q/in 
1 o/U/y4U 


Z 


pop2 


UU 1 4 


IVILAN/uK/zUUo/UU 1 4 


Crete-Heraklion 


LanL 


y© o/U/y4U 


J 


pop2 


UU 1 / 


IVILAIN/Un/zUU4/LiU 1 / 


Crete-Heraklion 


LanL 


no 07n/n/in 

yy o/u/y4(j 


5 


pop4 (sub4A) 


UU 1 o 


)VILAI\l/Un/zUU4/UIJ 1 o 


Crete-Heraklion 


LanL 


1 DZD 


7 O 

/y 


pop4 (sub4C^) 


UU 1 y 


IVILAN/uK/zUU4/UU 1 y 


Crete-Heraklion 


LanL 


1 DzD 


7 n 


pop4 (sub4C^) 




IVILAIN/uK/zUU4/LiUzU 


Crete-Heraklion 


LanL 


1 o/U/y4(J 


/ / 


pop2 




)VILAlN/UK/zUU4/UIJzo 


Crete-Lasithi/Sitia 


LanL 


1 Q/i n 
1 y4U 


z 


pop2 


uUz4 


KAr AK \ /ro /^r\r\A /rr\^A 
IVIL A n/ zUU4/ b Uz4 


Crete-Lasithi/Sitia 


LanL 


1 o/U/y4U 


z 


pop2 


KjUZd 


IVlLAN/uK/zUU4/bUzj 


Crete-Lasithi/Sitia 


LanL 


1 940 


2 


pop2 


kjUZo 


)VILAlN/Un/zUU4/UIJzD 


Crete-Lasithi/Sitia 


LanL 


1 Q7n/Q/in 
1 o/U/y4U 


z 


pop2 


uUo 1 


IVILAN/uK/zUUj/UUo I 


Crete-Lasithi/Sitia 


LanL 


1 o/U/y4U 


z 


pop2 


(jUoz 


IVILAN/uK/zUUj/LiUoz 


Crete-Lasithi/Sitia 


LanL 


1 940 


2 


pop2 




IVILAN/UK/zUUj/UUoo 


Crete-Lasithi/Sitia 


LanL 


1 o/U/y4U 


z 


pop2 




K/ir" A M /CD /ic\c\c, /r^n^c; 
IVILAlN/Lin/ zUUj/UUo J 


Crete-Lasithi/Sitia 


LanL 


QQ Q7n/Q/in 

yo o/U/y4U 


r 
J 


pop4 iSUD4Aj 




IVILAN/uK/zUUj/bUoD 


Crete-Lasithi/Sitia 


LanL 


QQ 07n/Q/in 

yy o/(j/y4(j 


5 


pop4 (sub4A) 


UU/ 1 


IVILAIn/UK/zUUd/UU/ I 


Crete-Lasithi 


LanL 


QQ Q7n/Q/in 

yo o/u/y4u 


5 
J 


pop2 


UU 1 Uz 


IVILAN/bn/zUUo/bU 1 Uz 


A -1-1- 1 / D -\ -Tl -\ ^ 

Attica (KaTinaj 


LanL 


QQ Q7n/Q/in 

y© o/u/y4u 


0 

J 


pop2 


GDI 04 


MCAN/GR/2006/GDl 04 


Crete-Heraklion 


CanL 


98 870/940 


20 


pop4 (sub4A) 


GD112 


MCAN/GR/2006/GD112 


Crete-Heraklion 


CanL 


98 870/940 


21 


pop4 (sub4A) 


GDI 54 


MCAN/GR/2008/GD154 


Central Macedonia-Thessaloniki 


CanL 


1 940 


22 


pop2 


GD167 


MCAN/GR/2008/GD167 


Thessaly-Trikala 


CanL 


98 870/940 


// 


pop2 


GDI 92 


MCAN/GR/2008/GD192 


South Aegean-Dodecanese/Rhodes 


CanL 


1 626/870/940 


22 


pop2 


GD201 


MCAN/GR/2008/GD201 


East Macedonia-Drama 


CanL 


1 870/940 


23 


pop2 


GD204 


MCAN/GR/2009/GD204 


East Macedonia-Kavala 


CanL 


1 870/940 


24 


pop2 


GD210 


MCAN/GR/2009/GD210 


Epirus-loannina/Dodoni 


CanL 


1 626 


25 


pop2 
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Table 1 Characteristics of the ^28 Leishmania infantum strains studied (Continued) 



uUz 1 1 


IVILAlN/bn/zUUy/uUz 1 1 


Peloponnese-Argolis 


LanL y 


o o/U/y4U 


5 
J 


pop2 


GD220 


IVlLAIN/uK/zUUy/UUzzU 


West Greece-Achaea 


LanL 


940 


26 


pop2 


KjUZZ I 


IVILAlN/Un/zUUy/UIJzz 1 


Epirus-Arta 


LanL 


DZD 


ZD 


pop2 




IVILAlN/UK/zUUy/UUzzo 


Ionian Islands-Corfu(Kerkyra) 


LanL 


OzD/y4U 


07 

z/ 


pop4 (sub4C^) 


GD229 


IVILAN/uK/zUUy/LiUzzy 


South Aegean-Cyclades 


LanL 


940 


/3 


pop2 




IVILAI\l/UK/zU 1 U/UUzo/ 


Central Macedonia-lmathia 


LanL 


OA n 
y4U 


OQ 
ZO 


pop2 


Anz 


IVInUIVI/AL/zUUD/Anz 


unknown 


\/i 1 

VL 


o/U/y4U 


JZ 


pop2 


A U/l 

An4 


IVInUIVI/AL/zUUD/An4 


unknown 


\/i 1 

VL 


o/U/y4U 


/ / 


pop2 


LtlVoyoz 


IVInUIVI/AL/zUUU/ J 1 oLLZn 


Ti rana 


\/l 1 
VL 


o/U/y4U 


JO 


pop4 (sub4B) 


AU 1 0 


\/ir"AM/Ai /onn/i /A Pii 
IVILAN/ AL/ zUU4/ AU 1 0 


unknown 


LanL 


n.Q. 


54 


pop2 


Lt)Vooy4 


IVILAI\l/AL/yo/L/oLj/4 


Ti rana 


LanL 


Q7n /Q/i n 
o/U/y4U 


c// 
-)4 


pop2 


1 1 1 
LU 


IVILAiN/ n n/zUUy/ LU 


Unknown 


LanL 


DZD 


4_) 


popo iSUDoD ) 


i\ie 


IVILAN/ H K/ zUUy/ IN 6 


Cetina(Brac) 


LanL 


DzD 


4D 


popo (SUDOD j 


I Vie 


IVILAIN/ n n/zUUy/ IVie 


Cetina(Makarska) 


LanL 


DzD 


4/ 


popo iSUDOD ) 


ivir 


IVILAIN/ HK/zUDo/IVir 


unknown 


LanL 


DzD 


4o 


pop4 (sub4B) 


Lib 


IVILAIN/ H K/ zUUy/Lio 


unknown 


LanL 


DzD 


/I n 

4y 


pop4 (sub4B) 


L 


IVILAIN/ riK/zUUy/L 


unknown 


LanL 


626 


47 


popo (SUDOD j 


1 1 
U 


IVILAIN/ n K/ zUUy/ U 


Dubrava(Korcula) 


LanL 


DzD 


4o 


pop4 (sub4B) 


ra 


IVILAIN/ H K/ zUUy/ r a 


Dubrava(Korcula) 


LanL 


DzD 


jU 


pop4 (sub4B) 


oU 


A M /UD/om 1 /cr^ 
IVILAIN/ n n/zU 1 l/oU 


Dubrovnik 


LanL 


DzD 


4_) 


popo iSUDoD ) 


A 1 

AL 


Kyir"AM/UD/oni i /ai 
IVILAIN/ n K/zU 1 l/AL 


opilt 


LanL 


DzD 


1 


popo iSUDOD ) 


Lb 1 


\/iur^\/i /CD /I nn^^/i cn/itc: 
IVIHUIVI/rK/ 1 yyo/LLIVI/ J 


Cevennes 


\/i 1 

VL 


DZD 


Dy 


popo iSUDoD ) 


LUz 


^/lu^^^yl /CD /I QQc: /I dmii/i 
IVInUIVI/rK/ lyyj/LrIN 1 14 


Cote d'Azur 


\/l 1 
VL 


DzD 


7n 
/U 


popo iSUDOD ) 


LU4 


^/lu^^^yl /CD /I QQ7 /I ci lo 
IVInUIVI/rK/ 1 yy//LoLzy 


Cevennes 


LL 


DzD 


7 1 
/ 1 


popo iSUDOD ) 


LLIVI- 1 1 oy 


\/ir" A M /cc /onm /i i \/i i i Qn 
IVILAIN/ Lo/zUUz/ LLIVI- 1 1 oy 


iDiza 


LanL 


DZD 


55 


popo iSUDoAj 


1 1 N /I 1 OOO* 

LLIVI- 1 zzo 


N /ir" A M /cc /onnQ /I i ^/l i ooo 
IVILAIN/ Lb/zUUo/ LLIVI- 1 zzo 


iDiza 


LanL 


626 


56 


popo (SUDoAj 


1 1 N/l 1 1A 1 * 

LLIVI- 1 z4 1 


N/ir" A M /CC /OnOQ /I 1 N/l 10/11 

IVILAIN/ Lo/zUUj/ LLIVI- 1 z4 1 


iDiza 


LanL 


DzD 


j/ 


popo iSUDoAj 


rlVII 


\/iur^\/i/cc/i nnQ /d\/i i 
IVIHUIVI/LV 1 yyo/rivi 1 


Majorca 


\/i 1 

VL 


DzD 


c:q 

JO 


popo iSUDoAj 


1 1 ^ /I no 1 * 
LLIVI-yo 1 


N/iur^N/i /cc /onni /i i ka ooi 
IVIriUIVI/Lb/zUU l/LLIVI-yo 1 


Majorca 


\/l 1 
VL 


626 


59 


popo (SUDoAj 


LLIVI- 1 UUo 


N/ir" A M /cc /onm /i i n/i inno 
IVILAIN/ Lo/zUU 1 / LLIVI- 1 UUo 


Majorca 


LanL 


DzD 


DU 


popo iSUDoAj 


LLIVI- 1 1 jU 


K/iuPiK/i /CQ /mm /I IK/1 1 1 c^n 
IVInUIVI/Lo/zUUz/ LLIVI- 1 1 jU 


Majorca 


r"i /'ui\/_l'\ 1 
LLinlV-l-j 


DzD 


/^i 
D 1 


popo iSUDoAj 


1 1 N /I 1 HQ 

LLIVI- 1 Uoo 


N/ir"A M /cc/onni /i i n/i ihqo 
IVILAIN/ Lb/zUU 1 /LLIVI- 1 Uoo 


Majorca 


LanL 


626 


62 


popo (SUDoD j 


LLM-1 116* 


MCAN/ES/2002/LLM-l 116 


Madrid 


CanL 1 


626 


63 


pop3 (sub3B^) 


LLM-1136* 


MCAN/ES/2001 /LLM-1 136 


Madrid 


CanL 1 


626 


64 


pop3 (sub3B^) 


LLM-1 181* 


MHOM/ES/2002/LLM-1181 


Madrid 


VL(HIV+) 1 


626 


65 


pop3 (sub3B^) 


IMT184* 


MHOM/PT/1993/IMT184 


Lisbon 


VL 1 


626 


65 


pop3 (sub3B^) 


IMT339* 


MCAN/PT/2003/IMT339 


Lisbon 


CanL 1 


626 


66 


pop3 (sub3B^) 


IMT205* 


MCAN/PT/1995/IMT205 


Alentejo 


CanL 1 


626 


65 


pop3 (sub3B^) 


IMT288* 


MHOM/PT/2002/IMT288 


Alto Douro 


VL 1 


626 


67 


pop3 (sub3B^) 


IMT316* 


MCAN/PT/2003/IMT316 


Lisbon 


CanL 1 


626 


68 


pop3 (sub3B^) 



reference strains included for comparison, the number of strains per geographical source is given in parenthesis; ^strains analyzed in previous studies 
[14,15,36,52,55,61,68]; TR, Turkey; CY, Cyprus; BG, Bulgaria; GR, Greece; AL, Albania; HR, Croatia; FR, France; ES, Spain; PT, Portugal; ^VL, visceral leishmaniasis; CanL, 
canine leishmaniasis; CL (HIV+), cutaneous leishmaniasis (co-infection with HIV); "^based on MLEE and K26 typing; n.d. not done; ^The 71 microsatellite profiles 
shared between different strains are indicated in italics; ^Population assignment at K = 4 by STRUCTURE analysis, sub-structuring upon re-analysis is shown in par- 
enthesis; infantum MON-1 clone isolated from a case of mixed infection with MON-1 and IVlON-37 strains [34,36]. 
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or eluted in TE buffer and stored at -20°C, whereas work- 
ing solutions were kept at 4°C until use. 

K26-PCR, MLEE and MLMT analysis 

The /<26-PCR assay was applied to discriminate MON-1 
from other L infantum zymodemes according to the length 
polymorphism of the K26 gene, as described previously 
[14]. MLEE analysis based on starch gel electrophoresis of 
fifteen enzymatic loci [6] was performed for all but six 
Turkish strains, at the University of Montpellier (France), 
to complement K26 typing. 

Multilocus microsatellite typing (MLMT) was applied 
using fourteen microsatellite markers previously found 
highly discriminative within the L. donovani complex 
[52,61]. Microsatellite fragments were amplified using 
established PGR protocols [52]. Amplicons were resolved 
on high-resolution MetaPhor® gels (3.5-4.5%) using a large 
gel horizontal electrophoresis system with built-in recircu- 
lation (Owl Large Gel System, Thermo Scientific) and 
screened for length polymorphisms using the Alphalma- 
ger® imaging system (Alpha Innotech Gorp.) as previously 
described [36]. Alternatively, PGR assays were performed 
using fluorescence conjugated primers (WellRed dyes, 
Proligo) and analyzed by capillary electrophoresis (GEQ™ 
8000, Beckman Goultier), as described elsewhere [36]. 
Gontrol strains were included in both PGR protocols and 
reading methods to confirm the reproducibility of the 
results. 

Population structure and phylogenetic analysis 

The MLMT dataset of 128 strains, including that of 30 
previously analyzed strains [15,36,52,55,61], was processed 
by models based on Bayesian statistics and genetic dis- 
tances. Model-based analysis was performed by STRUG- 
TURE v2.3.1 [62] to infer the population structure and 
identify genetically distinct L, infantum clusters. Based on 
multilocus genotypes and independent of a mutation 
model, this algorithm was used to determine patterns of 
allele frequencies per analyzed locus and estimated the 
most probable degree of differentiation. Individuals were 
assigned in discrete genetic populations and exposed their 
membership to each estimated population by genotypic 
fractions. The allele frequencies among populations were 
correlated by admixture modeling for a series of runs 
using a 'burn-in period of 2x10^ iterations and probability 
estimates based on 2x10^ of MGMG repeats. For each 
possible number of populations between K = 1 and K = 10, 
ten independent simulations were conducted to estimate 
the delta (A)K values [63]. The likelihood of population 
number was calculated by STRUGTURE HARVESTER 
[64], implementing the Evanno method [63]. Subsequent 
STRUGTURE analyses were performed to achieve assign- 
ment at different hierarchies. 



Phylogenetic analysis was performed based on microsat- 
ellite genetic distances, which were calculated by MIGRO- 
SAT software using the proportion of shared alleles 
distance measure (Das)- The resulting distance matrix 
input was processed by PHYLIP v3.6 to construct a 
Neighbour- joining (NJ) tree with confidence intervals by 
bootstrapping (1000 replications). A consensus tree with 
the corresponding bootstrap values and branch lengths 
was re-constructed using GENEIOUS v5.4 [65] and the 
derived mid-point rooted NJ tree was edited by FIGTREE 
(http://tree.bio.ed.ac.uk/software/figtree). Genetic related- 
ness and recombination were evaluated through Neigh- 
borNet networks constructed by SplitsTree4 software and 
through the population membership coefficients esti- 
mated by STRUGTURE. The MLMT data were analyzed 
by GDA (http://hydrodictyon.eeb.uconn.edu/people/plewis/ 
software.php) with respect to the proportion of poly- 
morphic loci (P), allelic variation (A), inbreeding coefficient 
(F/5), expected (HJ and observed (H^) heterozygosity. GDA 
analysis was performed per microsatellite locus and per 
population/subpopulation defined by STRUGTURE. F-sta- 
tistics were estimated for all populations/subpopulations to 
evaluate genetic differentiation and gene flow. For this, the 
F^T^ values (pairwise Wrights fixation index) with their cor- 
responding /7-values (confidence test) were calculated using 
MSA [66]. Moreover, the total strain-set was sorted by geo- 
graphical source to calculate the descriptive statistics per 
endemic country. Similarly, a separate descriptive analysis 
was performed for Turkish/Greek MON-1 and MON-98 
strains based on their zymodeme type and origin. 

The spatial distribution of the observed genetic struc- 
ture at population and individual level was assessed based 
on allelic similarity by factorial correspondence analysis 
(EGA) using GENETIX v4.03 [67]. 

Ethics statement 

The Leishmania strains designated with an asterisk (*, 
Table 1) and those originating from SW Europe (^, 
Table 1), which were used in this microsatellite analysis 
for comparison, have been the object of previous publica- 
tions [14,15,36,52,55,61,68]. 

Strains from Albania were obtained from the National 
Reference Gentre of Leishmania at MontpeUier, France. 
Ethical considerations are described in Petrela et al [45]. 
Strains from Greece were obtained from the already- 
existing ''Leishmania cryobank collection" of the Medical 
School of the University of Grete. Informed written con- 
sent was obtained from aU patients. Both the study and 
protocols used were approved by the Ethical Scientific 
Gommittee of the Medical School (University of Grete, 
Greece). Ethical considerations are described in Kuhls et al 
[14,15,36,52,55,61,68]. Strains from Turkey were collected 
in the context of previous studies, as indicated above. 
Human strains were obtained from an already-existing 
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"Leishmania Group cryobank collection". Informed written 
consent was obtained from all patients. The study was ap- 
proved by Local Animal Care and Ethics Committee of the 
School of Medicine Izmir, Turkey. 

All samples were coded and anonymized, and all Leish- 
mania strains were isolated as part of normal diagnosis. 
Dog samples were collected by veterinarians after the 
owners consent, following protocols that were approved 
by the Ethics Committees of the respective Institutions in 
Turkey, Cyprus, Greece and Croatia. Questionnaires with 
personal, epidemiological, and clinical data for each dog 
were also completed. 

Approval from the respective Institutions to use the 
samples for research purposes was received in all cases. 
The present study was also reviewed and approved by 
the Ethics Committee of the Hellenic Pasteur Institute, 
Athens, Greece. 

Results 

K26-PCR and MLEE analysis 

Of the total 128 strain-set, 125 strains were character- 
ized at species level and zymodeme type by /<r26-PCR, 
and 122 of them also by MLEE, in the context of either 
this or previous studies [15,36] (Table 1). A single 626 bp 
K26 amplicon characteristic for L, infantum MON-1 was 
presented by a few Turkish and Greek strains as well as by 
all strains from Cyprus and Croatia, including two clones 
of strain CD44 (CD44cl.4 and CD44cl.5) isolated from a 
canine with a mixed L, infantum MON-1 /L. donovani 
MON-37 infection in Cyprus [34,36]. The 940 bp/<26 
amplicon, which could correspond to either L infantum 
MON-1 or MON-98 zymodeme and was observed so far 
only in Greek strains, was presented by all strains from 
Turkey, Bulgaria and Albania, except for two Turkish 
strains from regions other than the Aegean. These were 
strains EPS 8 and EP134 isolated at the Black Sea and the 
Mediterranean region close to Cyprus, respectively. 
Besides the main 940 hp K26 amplicon observed for 
most strains from Greece, two Greek strains (GD192 
and GD223) gave multiple /<26-PCR products (626/870/ 
940 bp and 626/940 bp, respectively). In addition to 
seven previously typed L. infantum MON-98 strains 
from Crete, MLEE analysis identified 14 new strains of 
MON-98 from eight human and canine isolates from 
mainland regions of Greece and the island of Crete as 
well as six canine isolates from the Aegean region of 
Turkey (Ku^adasi) (Figure 1). Overall, our sample-set 
comprises 101 MON-1 strains, 21 MON-98 strains from 
Turkey and Greece, and 6 Turkish L infantum strains 
of unknown zymodeme type. Interestingly, 21 out of 70 
Turkish and Greek strains (30%) were identified as 
MON-98, which is the highest percentage of this zymo- 
deme found so far in Turkey and Greece. 



Analysis of microsatellite loci and correlation with 
zymodeme type 

The inter- and intra-zymodeme genetic diversity among 
1281. infantum MON-1 and MON-98 strains from south 
Europe was assessed by a panel of previously described 
microsatellite markers located on nine different chromo- 
somes [15,36,52,53,61]. Table 2 summarizes the character- 
istics and descriptive statistics per analyzed locus for the 
overall strain-set. 

Different degrees of polymorphism were observed 
among the analyzed loci. As shown in Table 2, the high- 
est allelic richness, expected and observed heterozygos- 
ities were demonstrated for markers Lm4TA, Lm2TG 
and Li22-35 (E), with the expected heterozygosity (mean 
He = 0.29) always being much higher than the observed 
(mean = 0.02). Inbreeding coefficients (F/5) were sub- 
stantially high (mean F/5 = 0.92), as four markers pre- 
sented absolute heterozygote deficiency (F/5 = 1), and 
marker Li46-67 (C) was monomorphic as previously ob- 
served [52,54]. 

Separate analysis was performed to investigate intra- 
zymodeme diversity and avoid possibly biased results 
due to the larger number of MON-1 strains in the total 
sample-set. The measures of genetic variability of Turk- 
ish/Greek strains demonstrated extensive inbreeding in 
both zymodeme populations and just trivial differences 
between MON-1 and MON-98 zymodemes (Table 3). 
However, some alleles at markers Lm2TG (22 repeats), 
CS20 (22 repeats), Li45-24 (G) (18 repeats) and LIST7039 
(16 repeats) were only found in MON-98 strains. 

Population structure inference 

The Bayesian model-based clustering algorithm imple- 
mented in the STRUCTURE software was run by in- 
creasing K values from K = 2 to K = 10 and with seven 
geographically-defined strain-groups, as strains from 
France (FR), Spain (ES) and Portugal (PT) (Table 1) repre- 
senting SW Europe were considered to constitute one 
population, according to earlier reports [15,36,52,55,61]. 

Calculation of AK indicated the existence of four dis- 
tinct genetic populations (Figure 2A). Population 1 (23 
strains, all MON-1) comprised all strains (canine iso- 
lates) from Cyprus (CY). Population 2 (51 strains) was 
formed by half of the strains from Turkey (8 TR), the 
single isolate from Bulgaria (BG) and the majority of 
strains from Greece (GR) and Albania (AL) (38 GR and 
4 AL). Five of these strains belonged to zymodeme 
MON-98 and were isolated from dogs in remote places 
of Greece, while all other were MON-1 strains. Popula- 
tion 3 (26 strains, all MON-1) incorporated all strains 
from SW Europe (19 FR/ES/PT), more than half of the 
strains from Croatia (6 HR) as well as the only TR strain 
isolated in the Black Sea region. Population 4 (28 
strains, MON-1 and MON-98) enclosed the remaining 
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Figure 1 Geographical distribution of the 128L infantum strains studied. Strains are represented by coloured circles that correspond to 
their population assignment at K = 4 by STRUCTURE analysis (Figure 2A); Population 1 -green, population 2-yellow, population 3-blue and 
population 4-red circles. The circles representing the Turkish strain isolated from Black Sea region ((^orum prefecture) and all strains from Croatia 
are colored according to their subpopulation assignment as in Figure 4. L. infantum MON-98 strains are indicated with an asterisk. Strains included 
for comparison from France, Spain and Portugal (Table 1) are not shown. 



TR (7), GR (16), AL (1) and HR (4) strains. Sixteen out 
of the 21 MON-98 strains from Ku^adasi (Turkey) and 
Crete (Greece) were assigned to this population. Popula- 
tion assignment per analyzed strain at K = 4 is desig- 
nated also in Figure 1 and Table 1. 

Population structure at K = 4 was well-established in 
all ten iterations and different genetic distance analyses 
(Figures 3, 4 and Additional file 1: Figure SI) corrobo- 
rated the existence of the four populations determined 
by the model-based analysis. However, the assignment of 
HR strains to POP3 and POP4 by STRUCTURE analysis 
at K = 4 (Figure 2A) was found inconsistent with their 
placement at the mid-point rooted NJ tree (Additional 
file 1: Figure SI) as well as with the distribution of the re- 
spective splits in the phylogenetic network (NeighborNet) 
(Figure 3). The FCA plot also confirmed the existence of 
at least four main populations (Figure 4). In overall, the 
three-dimensional representation mirrored the genetic 
isolation of CY strains (POPl), the differentiation of SW 



European strains (POPS) from all other populations as 
well as between the Balearic Islands and FR/mainland ES/ 
PT, and the diversification of TR/HR strains from SW 
European strains. 

The population structure defined by STRUCTURE 
was however, unstable for other K values and subdivision 
was detected by re-analyzing these four populations in- 
dependently. Population 1 split further in two subpop- 
ulations (not shown), but subdivision was not supported 
by other genetic distance analyses (Figures 3, 4 and 
Additional file 1: Figure SI). 

Population 2 was divided in three subpopulations 
(not shown); however there was no correlation between 
subpopulations and host specificity, as both human and 
canine isolates were found in every subcluster. Addition- 
ally, the assignment of strains to these three subpopula- 
tions was not related to their geographical origin as 
strains from both mainland and coastal areas of TR and 
GR never divided, and the GR strains always grouped 
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Table 2 Characteristics and descriptive statistics by 
microsatellite locus for the 128L infantum strains 
analysed in this study 



Locus 


Chr 


Repeat array 


P 


A 


He 


Ho 


F/s 


Lm2TG 


1 


TG (22-27) 


1 


6 


0.69 


0.09 


0.88 


Lm4TA 


1 


TA (10-17) 


1 


8 


0.76 


0.10 


0.87 


Li41-56 (B) 


36 


CA (9-11) 


1 


3 


0.37 


0.00 


1.00 


Li71-7 (R) 


30 


CA (12-14) 


1 


3 


0.12 


0.02 


0.87 


Li22-35 (E) 


1 


CA (10-13,15,16,18) 


1 


7 


0.64 


0.05 


0.92 


Li23-41 (F) 


25 


GT (13,15-17) 


1 


4 


0.35 


0.02 


0.96 


Li45-24 (G) 


36 


CA (14-18) 


1 


5 


0.51 


0.02 


0.95 


CS20 


19 


TG (18-20) 


1 


3 


0.29 


0.01 


0.97 


LIST7031 


10 


CA (11,12) 


1 


2 


0.11 


0.01 


0.93 


LIST7039 


30 


CA (14-16) 


1 


3 


0.09 


0.01 


0.92 


Li46-67 (C) 


31 


CA (9) 


0 


1 


0.00 


0.00 


0.00 


Li71-33 (P) 


31 


TG (10-12) 


1 


3 


0.03 


0.00 


1.00 


Li71-5/2 (Q) 


35 


CA (8-10) 


1 


3 


0.08 


0.00 


1.00 


Tube A 


34 


CA (9) 


0 


2 


0.01 


0.00 


1.00 


Overall 






0.86 


3.79 


0.29 


0.02 


0.92 



Locus, analysed microsatellite marker; Chr, chromosome location; Repeat array 
per locus, in parenthesis are the overall dinucleotide repeat numbers 
identified; P, proportion of polymorphic loci; A, number of alleles per locus; Hg, 
Nei's unbiased expected heterozygosity; Ho, observed heterozygosity; F/s, 
inbreeding coefficient. 



with those from BG and AL. Moreover, no differenti- 
ation based on zymodeme type was detected between 
GR MON-98 strains and BG/GR/AL MON-1 strains. In 
addition, the subdivision of POP2 was not supported by 
the different distance-based analyses used (Figures 3, 4 
and Additional file 1: Figure SI). 

Population 3 was divided in three stable subpopula- 
tions (Figure 2B), All but one strain from the islands of 
ES were assigned to a distinct subpopulation (sub3A), as 
previously observed [15]. Sub3B^ enclosed all strains 



Table 3 Genetic diversity between and within MON-1 and 
MON-98 strains from Greece and Turkey based on the 
MLMT profiles in 14 analyzed markers 



Population 


n 


Matches 


P 


A 


He 


Ho 


F/s 


TR/GR MON-1 


44 


2 


0.71 


2.64 


0.19 


0.02 


0.87 


GR MON-1 


40 


3(5) 


0.57 


2.14 


0.17 


0.02 


0.86 


TR MON-1 


4 


0 


0.64 


1.93 


0.24 


0.03 


0.88 


TR/GR MON-98 


21 


0 


0.57 


2.21 


0.15 


0.03 


0.82 


GR MON-98 


15 


1 


0.43 


1.93 


0.17 


0.02 


0.89 


TR MON-98 


6 


1 


0.29 


1.29 


0.05 


0.05 


0.00 



Nine identical genotypes (matches) were identified between MON-1 and 
MON-98 strains from Turkey (TR) and Greece (GR). The number in parenthesis 
includes two common genetic profiles found between MON-1 and MON-98 
strains from Greece, n, number of analyzed strains per population; P, propor- 
tion of polymorphic loci; A, number of alleles per locus; Hg, Nei's unbiased ex- 
pected heterozygosity; Ho, observed heterozygosity; F/s, inbreeding coefficient. 



from FR, PT and the mainland of ES as well as single 
strains from TR, HR and the ES island of Majorca. Of 
note is that the three strains from TR, HR and Majorca 
that grouped with strains from the mainland ES in 
sub3B^ (Pi gure 2B), were placed in intermediate posi- 
tions at the Splitstree and PCA plots (Pigures 3 and 4) 
and that the single strains EPS 8 and AL from TR and HR 
respectively, formed a separate subgroup distant from 
other strains of POPS at the NJ tree (Additional file 1: 
Pigure SI). Moreover, strains EP88 and AL presented large 
splits in the phylogenetic network (Pigure 3) and were 
found more closely related to a group of TR/GR isolates 
(strains EP134, GD18 and GD19) that were assigned at 
POP4 by STRUCTURE. This finding is also illustrated in 
the PCA plot (Pigure 4). Accordingly in the NJ tree, strains 
EP88 and AL branched out and clustered with the same 
group of TR/GR strains after the formation of POPl 
(Additional file 1: Pigure SI). All other strains from HR 
(Lu, Ne, Me, SO and C) grouped in a separate homoge- 
neous subpopulation (sub3B^) (Pigure 2B), which was 
supported by all analysis methods (Pigure 3, 4 and 
Additional file 1: Pigure SI). 

Re-analysis of population 4, which included the majority 
of the MON-98 strains, revealed further sub-structuring 
(4 stable subpopulations) based on both zymodeme type 
and geographical origin (Pigure 2C). Sub4A included only 
MON-98 strains from Turkey and Greece, while the Split- 
stree plot (Pigure 3) illustrated gene flow between them. 
The group of HR strains (Pa, Mr, GS and U) and one strain 
from Albania (LEM3962), which was assigned to sub4B, 
presented a reticulate pattern in the Splitstree network 
(Pigure 3) rather closely related to two TR strains (EP58 
and EP68) which were assigned to POP2. This was also 
apparent in the NJ tree (Additional file 1: Pigure SI) and 
the PCA plot (Pigure 4). Turkish and Greek MON-1 and 
MON-98 strains were found in sub4C^, whereas sub4C^ 
enclosed only GR MON-1 strains, most of which pre- 
sented the 626 hp K26 amplicon (Table 1). 

Genotypic characterization and F-statistics for populations 
and subpopulations 

The measures of genetic diversity calculated for each main 
population defined by STRUCTURE are summarized in 
Table 4. Out of the 128 analyzed strains, 50 strains pre- 
sented unique microsatellite profiles and 78 strains shared 
the same MLM Type with at least another strain. Identical 
profiles were found only among strains that were assigned 
in the same main population (Tables 1 and 4). The genetic 
polymorphism exhibited per analyzed locus for each popu- 
lation is given in Additional file 2: Table SI. 

All strain populations were polymorphic at markers 
Lm2TG, Lm4TA, Li22-35 (E), Li23-41 (P) and CS20. The 
highest degree of polymorphism and allelic richness in all 
markers was observed within POPS and POP4, where 
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Figure 2 Population structure of the 128L infantum strains inferred by STRUCTURE analysis. Population structure is based on tine MLMT 
data observed for eacli analysed strain at 14 microsatellite markers. Each population is represented by different colour and differently coloured 
segments in the vertical bar correspond to the membership of each strain to different populations. The length of each coloured segment is 
proportional to the estimated membership coefficient of each strain per population. The percentage of strains enclosed per population/ 
subpopulation based on their geographical origin (TR, Turkey; CY, Cyprus; BG, Bulgaria; GR, Greece; AL, Albania; HR, Croatia; FR, France; ES, Spain; 
PT, Portugal) and L infantum zymodeme (MON-1, MON-98), is demonstrated by large and small pie-charts, respectively. A. The likelihood of 
population number was calculated by STRUCTURE HARVESTER [64] implementing the Evanno method [63]. The major peak at K = 4 at the AK 
graph indicates the existence of four populations (populations 1-4) for the total strain set, which were constantly observed in all independent 
simulations (10 runs). Population subdivision was investigated by re-analyzing separately each population at K = 4; subdivision of populations 1 
and 2 is not shown. B. Three geographically determined subpopulations (subs 3A, 3B^ and 3B^) were observed upon re-analysis of population 3. 
Bar plots for both K = 2 and K = 3 are shown. C. Four stable subpopulations (subs 4A, 4B, 4C^ and 4C^) based on both geographical origin and 
zymodeme type, were identified for population 4. Bar plots for both K = 3 and K = 4 are shown. 



markers LIST7039 and Li71-33 (P) were also polymorphic. 
Comparable variation was detected within POP2, al- 
though this population included twice the number of 
strains comprising POPS and POP4, In contrast, POPl, 
exclusively formed by CY strains, presented the lowest 
polymorphism and allelic variation. Strains within POP2 
and POP3 further segregated at markers LIST7031 and 
Li71-5/2 (Q). Similarly, marker Li41-56 (B) was poly- 
morphic only among strains in POP2 and POP4. All main 
populations were characterized by much lower observed 
heterozygosities than the expected ones (mean = 0.17 
and mean = 0.02) as well as high inbreeding coeffi- 
cients (mean ¥js = 0.86) (Table 4). 

Regarding the subpopulations defined by STRUC- 
TURE and largely supported by the genetic distance 



analyses, the lowest levels of polymorphism, heterozy- 
gosity and inbreeding were found in sub3B^ (HR strains) 
and sub4B (AL/HR strains) (Additional file 3: Table S2). 
The fact that lower F/5 values are observed at subpopu- 
lation level indicate that Wahlund effect (population 
subdivision) accounts at least partially to the high in- 
breeding coefficients found for the main populations. 

The four populations identified by model- and distance- 
based analyses were further supported by F-statistics 
(Table 5) indicating strong differentiation between them. 
CY strains that grouped in POPl, demonstrated the great- 
est genetic isolation from all other populations (F57^^0.546- 
0.655). As shown in Additional file 4: Table S3, the 
observed subpopulations were also strongly supported by 
F-statistics (p < 0.006). 
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TR/GRMON-98 



Figure 3 Phylogenetic network (NeighborNet) constructed by Splitstree for the 128L infantum strains studied. The NeighborNet 
algorithm detected the same population substructures inferred by STRUCTURE (Figure 2, Table 1) for population 3 (Sub3A, Sub3B^ and Sub3B^) 
and population 4 (Sub4A, Sub4B, Sub4C^ and Sub4C^). Subdivision of populations 1 and 2 is not supported by NeighborNet analysis. TR, Turkey; 
CY, Cyprus; BG, Bulgaria; GR, Greece; AL, Albania; HR, Croatia; FR, France; ES, Spain; PT, Portugal; MON-1, L infantum zymodeme MON-1; MON-98, L 
infontunn zymodeme MON-98. Strains that presented identical MLM Types (Table 1) are underlined. Strains are represented by coloured boxes that 
correspond to their population or subpopulation assignment according to STRUCTURE analysis (Figure 2). 



Correlation of MLMT profiles with origin, zymodeme type 
and host 

The strain-set was also grouped based on geographical 
source and their MLM Types were analyzed to compare 
the degree of polymorphism between and within each 
endemic country. At least one matching pair was de- 
tected within each endemic country (Table 1). 

Notably, the same MLMT profile (No 29) was identi- 
fied for 17 out of 23 canine isolates from different pre- 
fectures of Cyprus and a second one (No 33) was shared 
by two canines from Paphos prefecture. Seven distinct 
profiles (No 45-51) were identified for the ten analyzed 
strains from Croatia. Of these, three multilocus geno- 
types (No 45, 47 and 48) were observed in more than 
one Croatian strain. Two canine isolates from Albania 
were indistinguishable by MLMT (No 54). The same 
MLMT profile (No 37) was identified among MON-98 
canine isolates from Ku^adasi district (Aegean region), 
whereas a second one (No 41) was shared by strains of 



unknown zymodeme (EP58 and EP68) that were isolated 
from dogs in coastal areas of the Aegean (Mugla and 
Ku^adasi, respectively). 

Only 20 out of 54 strains analyzed from Greece pre- 
sented distinct MLM Types (Table 1). The same MLMT 
profile (No 2) was shared by MON-1 strains (human 
and canine isolates) from the capital of Athens, the is- 
land of Mytilene and Lasithi prefecture in Crete. An 
identical profile (No 22) was shared by MON-1 strains 
(canine isolates) from Rhodes island and the capital of 
North Greece Thessaloniki, where the second largest 
seaport is located. Another two MLMT profiles (No 5 
and 19) were identical between either MON-98 strains 
or MON-1 strains from Crete. MLM Type No 3 was 
common between MON-1 strains (human isolates) from 
Athens and MON-98 strains (canine isolates) from a 
harbor area in Attica, northern Peloponnese and Crete. 

Furthermore, identical MLM Types were identified be- 
tween strains originating from different countries, except 
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Figure 4 Spacial distribution of the 128L infantum strains by factorial correspondence analysis (FCA). Subdivision of population 1 (green 
circle) and population 2 (yellow circle) is also not demonstrated in the FCA plot. The subpopulations observed at STRUCTURE (Figure 2B and 2C) 
and NeighborNet analysis (Figure 3) for population 3 (blue circle) and population 4 (red circle) correspond to those illustrated by FCAnalysis. 
These are Sub3A (dark blue squares), Sub3B^ (light blue squares), Sub3B^ (blue squares), Sub4A (dark pink squares), Sub4B (light pink squares), 
Sub4C^ (dark red squares) and Sub4C^ (red squares). TR, Turkey; CY, Cyprus; BG, Bulgaria; GR, Greece; AL, Albania; HR, Croatia; FR, France; ES, Spain; 
PT, Portugal; MON-1, L infantum zymodeme MON-1; MON-98, L infantunn zymodeme MON-98. Strains sharing the same MLM Type (Table 1) are 
not shown. The colours of the circles designating each subpopulation correspond to those in Figures 2B, 2C and Figure 3. 



Table 4 Descriptive statistics by population (POP) and 
comparison of the F/s values between the four main 
populations defined by STRUCTURE analysis 



POP 


Origin 


Total no of 
MLMT profiles 


P 


A 


He 


Ho 


F,s 


1 (23) 


CY 


6 


0.36 


1.36 


0.03 


0.01 


0.72 


2(51) 


TR, BG, GR, AL 


23 


0.64 


2.29 


0.12 


0.02 


0.87 


3(26) 


TR, HR, FR, 
ES, PT 


22 


0.86 


2.71 


0.28 


0.02 


0.91 


4(28) 


TR, GR, AL, HR 


20 


0.71 


2.86 


0.24 


0.04 


0.82 


Mean 






0.64 


2.30 


0.17 


0.02 


0.86 



The data are based on the analysis 128L infantum strains by 14 microsatellite 
markers. The number of strains assigned per population is shown in 
parenthesis. Origin acronyms indicate the strain geographical source. P, 
proportion of polymorphic loci; A, number of alleles; Hg, Nei's unbiased 
expected heterozygosity; Hq, observed heterozygosity; F/s, 
inbreeding coefficient. 



for strains from Croatia and Cyprus. Notably, the single 
canine isolate from Bulgaria shared the same MLMT 
profile (No 11) with one Albanian and three Greek 
strains. Also, two MLM Types (No 13 and No 25) were 
common between strains from Greece and Turkey. One 
of them was found in canine isolates from the Aegean 
region (Cyclades-Greece and Bodrum-Turkey, respect- 
ively) and a human isolate from Central Greece. The 
other MLM Type was shared solely by MON-1 canine 
isolates from north inland regions of Greece and Turkey 
and one strain of unknown zymodeme from Mugla 
(Aegean region). As shown in the Splitstree plot, these 
strains were separated in a clonal subgroup (Figure 3). 

Characteristic alleles only for CY isolates were identi- 
fied at markers Lm2TG and Lm4TA. Also, HR strains 
presented distinct alleles at marker Li71-7 (R) and GR 
strains at markers Li22-35 (E) and Li23-41 (F). No cor- 
relation was observed between specific microsatellite 
profiles and human or canine hosts. 
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Table 5 Fst values and corresponding p-values for the 
populations identified at K = 4 



Populations 


1 


2 


3 


4 


1 (23) 


0 


0.655 


0.546 


0.558 


2(51) 


<0.001 


0 


0.526 


0.409 


3(26) 


<0.001 


<0.001 


0 


0.358 


4(28) 


<0.001 


<0.001 


<0.001 


0 



Population assignment is based on STRUCTURE analysis; In parenthesis, the 
number of strains comprising each population; Fst, pairwise Wright's fixation 
index (upper right triangle); p-values, confidence test (lower left triangle). 



Discussion 

The aim of the present study was to evaluate epidemio- 
logical aspects of L infantum causing VL in the SE edge 
of Europe focusing on Turkey, Cyprus and Greece but 
also including strains from Balkan countries. 

Initially, the /<26-PCR assay was applied to confirm 
the presence of L. infantum spp. and simultaneously to 
distinguish MON-1 from non MON-1 zymodemes within 
the strain-set [14]. Application of the K26 typing tool, 
which discriminates between different L. infantum and 
L. donovani zymodemes, was also essential due to the 
recent detection of L donovani MON-37 in Cyprus [35] as 
well as the identification of the novel L donovani MON- 
308 and MON-309 zymodemes in Turkey [36]. Notably, 
strains belonging to these three L donovani zymodemes 
were found by MLMT to form a novel L donovani sensu 
lato group [36]. 

In consistence with previous findings [14], K26 typing 
revealed that all strains from SE Mediterranean and Bal- 
kan countries were L infantum presenting single K26 
amplicons of 626 bp or 940 bp, except for two Greek 
MON-1 strains showing both /<26-PCR products. In line 
with K26 typing, MLEE analysis confirmed the predom- 
inance of zymodeme MON-1. Of the 70 strains isolated 
from the Aegean region of Turkey, mainland regions of 
Greece and the island of Crete, 21 (30%) were typed, 
however, as zymodeme MON-98. This is the highest per- 
centage of MON-98 strains found so far [5,14,69] and it 
should be noticed that Greek MON-1 and MON-98 
strains were similarly distributed in the country. Since 
these are the only zymodemes found so far in the Balkans, 
the K26-PCR assay is adequate for Leishmania typing in 
this endemic region. 

The MLMT profiles obtained for 109 L infantum MON- 
1 and MON-98 strains from six SE European countries 
were compared to those of 19 previously analyzed L infan- 
tum MON-1 strains from SW Europe (France, Spain and 
Portugal) [15,52,61]. A small number of Greek MON-1 
and MON-98 strains and two Turkish L infantum strains 
of unknown zymodeme type (strains EP3 and EP16) 
(Table 1) had been also analyzed earlier [15,36,52,55,61]. In 
this study we have investigated additional strains from VL 
foci in Turkey and Greece and also included L infantum 



strains from Cyprus, Bulgaria, Albania and Croatia, for 
which genetic diversity and population structure had not 
been studied until now. By this we could confirm the exist- 
ence of substantial differentiation between L infantum 
strains from the eastern and western parts of south Europe, 
and that population structuring is geographically deter- 
mined, as previously demonstrated [70]. These findings are 
also in agreement with previous data according to which, 
Cypriot and Iberian strains could be distinguished upon 
processing their /<26-PCR products by RFLP [14]. 

Different types of population genetic analyses, including 
Bayesian inference (as implemented in STRUCTURE), 
distance-based (NJ and Neighbor Net in SplitsTree) and 
FCA revealed the existence of four putative populations. 
Population 1 was well-defined by all methods and con- 
sisted exclusively of the Cypriot strains of L infantum and 
was clearly separated from the other populations (Figures 3 
and 4). As expected, all SW European strains formed a 
distinct population (POP3). Bayesian statistics has also 
assigned one strain from the Black Sea area in Turkey and 
six strains from Croatia to this population that was, how- 
ever, not well-supported by the distance-based and FCA 
analyses, which indicates a rather intermediate position of 
these strains. Three subpopulations became apparent by 
separate re-analysis. The differentiation between strains 
from the Balearic Islands (sub3A) and strains from the 
Iberian Peninsula and France (sub3B^) confirms previous 
analysis [15,52,61]. All but one strain from Croatia that 
grouped in POP3, were assigned to the third subpopulation 
(sub3B^). The Turkish and Greek strains split among two 
populations, POP2 and POP4 with no correlation to their 
place and year of isolation, the zymodeme type, host speci- 
ficity or pathology. POP2 was composed by most strains 
from Turkey, Greece, Albania and the single strain from 
Bulgaria. Nevertheless, distance and FCA analyses revealed 
that this population was rather homogenous and was 
not divided into subpopulations. In contrast, POP4 that 
enclosed some Turkish/Greek, one Albanian and four 
Croatian strains, appeared to be highly diverse and con- 
sisting of four subpopulations. Sub4A was formed only 
by Turkish/Greek MON-98 strains and sub4B enclosed 
only Albanian/Croatian strains. The other Turkish/Greek 
strains, whether of zymodeme MON-1 or MON-98, were 
assigned to sub4C^ and sub4C^. 

In general, model- and distance- based algorithms iden- 
tified the same population and subpopulation divisions. 
The statistical support for the midpoint-rooted NJ tree 
was assessed by bootstrap analysis. However, sufficient 
bootstrap values were only observed for nodes between 
groups of strains found within the main clusters or sub- 
clusters (Additional file 1: Figure SI), while a spider-net 
pattern and large splits were often observed in Splitstree 
(Figure 3), probably due to the existence of homoplasy, 
intermediate or mixed genotypes and considerable genetic 
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migration. This was expected, since the analyzed strains 
belong to genetically close L. infantum zymodemes and 
were isolated from geographically close countries with fre- 
quent movements of host populations and could explain 
the discrepancies found in the assignment of some groups 
of strains or individual strains by different model- and 
distance-based approaches. Analogous congruity has been 
also observed elsewhere [15,52]. 

Compared to previous studies [15,36,52,55,61], we have 
now analyzed a much larger number of MON-1 and 
MON-98 strains from Turkey and Greece, but as before, 
we were unable to differentiate these strains according to 
their geographical origin, Turkey vs. Greece or Greek 
mainland vs. the island of Crete, or between MON-1 and 
MON-98 strains. Most MON-98 strains from both Turkey 
and Greece were assigned to the highly diverse POP4 
together with strains of zymodeme MON-1 from the 
same areas, but clustered in the two subpopulations 
sub4A and sub4C^ (Figure 2C). A distinct subpopula- 
tion of the MON-98 strains that grouped in POP2 was 
not supported by distance-based analyses. The reticulate 
network in Splitstree (Figure 3) clearly points to exten- 
sive gene flow and/or hybridization occurring among 
Turkish and Greek strains of both zymodemes. Our results 
were not affected when a group of representative non 
MON-1 strains, other than those of zymodeme MON-98, 
was included in our analysis (not shown). This finding is in 
agreement with previous analysis [15,36,52,55,61], accord- 
ing to which a small group of analyzed MON-98 strains 
from Crete did not cluster with other non MON-1 strains 
from SW Europe but always remained part of the MON-1 
population. Our observations reflect the high polymorph- 
ism observed among different groups of L. infantum 
MON-1 strains rather than indicate differentiation between 
MON-1 and MON-98 zymodemes. 

Identical MLM Types were found between MON-1 
strains from Turkey and Greece, which along with the ob- 
served shared membership of POP2 and POP4 and the sig- 
nificantly high inbreeding level within TR/GR MON-1 
(F/5 > 0.873) indicate considerable recombination and gene 
flow among them as previously suggested [15]. These find- 
ings may be associated with intensive host migration over 
long periods of time and common eco-epidemiological fea- 
tures, including a similar sand fly fauna [71], in Turkey and 
Greece. The detection of identical MLMT genotypes be- 
tween Greek strains of both zymodemes in harbor areas of 
the mainland and in Greek islands suggests parasite trans- 
fer between the mainland and the islands. 

The assignment of Greek human and canine isolates to 
the same genetic populations indicates the absence of host 
specificity. Since we have included a much larger number 
of strains, this observation enhances previous findings [15]. 

The genetic structure of strains from Crete appears to 
be a micrograph of the overaU epidemiology of VL in 



Greece, which is in contrast to the clear division be- 
tween strains from the mainland and islands of Spain, 
Majorca [8] and Ibiza [15]. At high hierarchy, the overaU 
high genetic diversity, significant gene flow and popula- 
tion subdivision have most probably resulted from suc- 
cessive recombination or even hybridization events. In 
paraUel, the observed low variability measures and con- 
siderable inbreeding observed within smaU subgroups of 
MON-1 or MON-98 strains from either inland regions 
of Greece or the island of Crete (Table 3) can be explained 
by successive subdivision followed by genetic drift within 
each of these groups and could indicate a founder effect. 
The existence of related clones of MON-1 as well as 
MON-98 strains could indicate selection of specific micro- 
satellite genotypes, as previously observed for MON-1 
strains from Madrid or Majorca and Ibiza [72]. 

L. infantum MON-1 canine strains from Cyprus were 
found genetically isolated from aU other genetic groups 
and formed a monophyletic group (Additional file 1: 
Figure SI) presenting the lowest degree of polymorphism, 
allelic deprivation and traces of heterozygosity, which in- 
dicate genetic drift. Seventeen out of the 23 strains shared 
the same MLMT genotype suggesting an epidemic out- 
break due to the same strain which has spread to different 
parts of the island. As expected the clone of strain CD44 
(CD44cl.5) shared common alleles and grouped with the 
other MON-1 strains from Cyprus. The identification of 
different allelic variants than those typically found for L. 
donovani MON-37 strains [35], responsible for both 
HVL and CL, suggests the absence of recombination or 
hybridization between MON-1 and MON-37 strains. 
The epidemiology of CanL in Cyprus has been doubtless 
affected by the eradication campaigns held against mal- 
aria between 1940 and 1950 which have massively re- 
duced the vector populations on the island. Another 
campaign was performed against echinococcosis be- 
tween 1970 and 1975 which diminished the dog popula- 
tions and altogether eliminated CanL for over 20 years 
[11]. Most probably this has caused severe bottleneck 
events in Leishmania populations leading to a limited 
parasite genetic pool characterized by significant inbreed- 
ing and genetic homogeneity. However, after the sand fly 
fauna has gradually recovered and CanL re-emerged, 
some degree of genetic variation was observed in the para- 
sites as a result of adaptation to selective pressure. 

The partition of canine MON-1 strains from Croatia 
between both SE (POP4) and SW (POP3) European pop- 
ulations (Figures 2, 4 and Additional file 1: Figure SI) 
mirrors the geographical position of the country. Despite 
the relatively small number of strains, we could document 
substantial genetic diversity for these strains. Croatia is sit- 
uated between East and West Europe and parasites can be 
introduced from its neighboring countries. The microcli- 
mate, the presence of vectors and reservoirs favor the 
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establishment and circulation of parasites in the Dalmatian 
coastal regions. The exportation from Croatia to neighbor- 
ing countries is equally possible, as adult VL cases acquired 
at the Dalmatian coast were registered in Austria and 
Switzerland [40,42]. The genetic differentiation observed 
between Croatian and SW European strains on one side 
and between Turkish, Greek and Croatian strains, on the 
other side, can be also seen by the distinct and intermedi- 
ate placement of the Croatian strains in the Splitstree 
(Figure 3). The rather low level of inbreeding and the 
substantially low variability values observed within these 
subpopulations (Tables 4 and Additional file 4: Table S3) 
indicate genetic drift, most probably caused by population 
substructure (Wahlund effect). Other explanations could 
be the restriction of the parasites' gene pool due to mass 
insecticide spraying for malaria control or the diversifica- 
tion of autochthonous cases circulating in Croatia for a 
long time [42]. Considering the unexpected genetic diver- 
sity of Croatian strains isolated from a rather small terri- 
tory, future investigations should use a larger sample set 
including human and sand fly isolates, and should include 
strains from the northern Adriatic region and Italy. 

Strains from Albania did not form a distinct genetic 
group in our analyses, but grouped either with strains 
from Turkey/Greece or Croatia. In contrast, a consider- 
able amount of genetic flow was detected obviously 
reflecting the frequent human migration and the unre- 
strained movement of canine and sand fly populations 
between these countries. HVL in Albania still remains 
an infantile disease and presents by far the highest mor- 
bidity rates in Europe [45] although relapses and resist- 
ance to drug treatment are only rarely reported [44-46] . 
It raises concerns that 57.1% of seropositive HVL cases 
hospitalized in the loannina prefecture (North Greece) 
are immigrants from the neighboring villages in south 
Albania and that the highest prevalence of clinically af- 
fected HVL (42.8%) and CanL (82.8%) in Greece is re- 
ported in the Thesprotia prefecture (Epirus) close to the 
border between Greece and Albania [31]. These findings 
are alarming for both countries and suggest that the dy- 
namics and complexity of VL in Albania could be simi- 
lar to those observed in Greece. CanL is endemic across 
the country with an estimated 16-17% seroprevalence, 
while sand fly species Phlebotomus neglectus, Phleboto- 
mus perfiliewi and Phlebotomus tobbi are widespread. It 
should not be underestimated that the maximum inci- 
dence of P, neglectuSy which is highly present in all Greek 
islands except Gavdos [5,71] and is a proven L, infantum 
vector in the island of Corfu neighboring the Albanian 
coast [44], is reported in endemic foci presenting the high- 
est VL morbidity rates. It should also be noticed that P, 
tobbi spp. despite its overall low abundance, is solely 
found in Albanian districts where VL morbidity exceeds 
1/1000 [44]. 



The lack of isolated Leishmania strains from Bulgaria 
has so far restricted the investigation of the population 
structure and dynamics. We have included in our study 
the first strain that was isolated from a dog in the 
Blagoevgrad province (southwestern Bulgaria) that bor- 
ders the Macedonia region of Greece (Figure 1). Our find- 
ings confirm the existence of at least L. infantum MON-1 
zymodeme in Bulgaria [48] and demonstrate that the 
MLM Type of this particular canine isolate is identical to 
that of three strains from Greece and one strain from 
Albania. This is not surprising as Bulgaria shares common 
geo-ecological characteristics with its neighboring coun- 
tries. The presence of L, infantum sand fly vectors and 
stable parasite circulation in Stara Zagora [49] highlights 
the necessity to conduct a comprehensive entomological 
and epizootic survey for the identification of Leish- 
mania vectors and the isolation of Leishmania parasites 
from sand flies, human and canine hosts. This would 
give the opportunity to thoroughly investigate the popu- 
lation structure of L, infantum spp. and disease dynam- 
ics in Bulgaria. 

Conclusions 

This is the first study that gives an insight into the popu- 
lation structure of L, infantum in Cyprus located at the 
SE edge of Europe and the Balkans (Turkey, Bulgaria, 
Greece, Albania and Croatia). L, infantum MON-1 was 
found to be predominant in this region, while 16.8% of 
strains were typed as MON-98. Hence, the K26 typing 
tool proves adequate for typing strains from SE Europe 
and the Balkans. Herewith, we have identified distinct 
genetic populations that reveal the clear differentiation 
between SE and SW European strains, the intermediate 
placement and further diversification of Croatian strains 
from SW strains and the formation of a monophyletic 
group enclosing all canine MON-1 isolates from Cyprus. 
The latter could be explained by the natural isolation of 
Cyprus being an island and a bottleneck effect caused by 
the subsequent eradication campaigns against malaria 
and echinococcosis that were held in Cyprus. Also, we 
have analysed the highest number of MON-98 strains 
found so far in Turkey and Greece and did not identify a 
distinct population including all MON-98 strains. The 
latter finding is strengthened by the separate analysis of 
MON-98 and MON-1 strains, which shows that the gen- 
etic distances observed between MON-98 and MON-1 
subgroups are similar to those found between different 
MON-1 subgroups. Altogether our results confirm and 
extend the study by Kuhls et al [15] using a larger strain- 
set from Turkey and Greece and including strains from 
other endemic Balkan countries. Notably, the clear diver- 
sification between Spanish strains from the mainland and 
the islands of Majorca and Ibiza is contrary to the cluster- 
ing of Greek strains from the mainland and the islands of 
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Corfu, Crete, Mytilene and Rhodes. Overall, the observed 
structuring is associated with the sand fly fauna, geo- 
ecological characteristics, historical background and im- 
portant events that have influenced the epidemiology of 
visceral leishmaniasis in these endemic countries. 

Additional files 



Additional file 1: Figure SI. Midpoint rooted Neighbor-joining tree 
constructed for the 128 L. infantum strains studied. The midpoint rooted 
tree is based on the Dps-distances calculated for the MLMT data at 14 
microsatellite markers for the total strain set. Bootstrap values only above 
50% are indicated at key nodes. L. infantum zymodemes, geographical 
origins and populations as inferred by STRUCTURE at K=4 (Figure 2A) are 
shown in colored boxes next to the tree. The colors designating each 
population correspond to those in Figures 2A, 3 and 4. TR, Turkey; CY, 
Cyprus; BG, Bulgaria; GR, Greece; AL, Albania; HR, Croatia; FR, France; ES, 
Spain; PT, Portugal. 

Additional file 2: Table SI. Descriptive statistics by microsatellite locus 
for each of the main four populations defined by STRUCTURE analysis. 

Additional file 3: Table S2. Descriptive statistics and comparison of F/s 
values between the sub-populations defined by STRUCTURE. 

Additional file 4: Table S3. Est values and corresponding p-values for 
the sub-populations identified by STRUCTURE re-analysis. 
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